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ABSTRACT 

Context. Observations of the y-ray sky with Fermi led to significant advances towards understanding blazars, the most extreme class 
of Active Galactic Nuclei. A large fraction of the population detected by Fermi is formed by BL Lacertae (BL Lac) objects, whose 
sample has always suffered from a severe redshift incompleteness due to the quasi-featureless optical spectra. 

Aims. Our goal is to provide a significant increase of the number of confirmed high-redshift BL Lac objects contained in the 2 LAC 
Fermi/LAT cataloge. 

Methods. For 103 Fermi/LAT blazars, photometric redshifts using spectral energy distribution fitting have been obtained. The pho- 
tometry includes 1 3 broad-band filters from the far ultraviolet to the near-IR observed with Swift/UVOT and the multi-channel imager 
GROND at the MPG/ESO 2.2m telescope. Data have been taken quasi-simultaneously and the remaining source-intrinsic variability 
has been corrected for. 

Results. We release the UV-to-near-IR 13-band photometry for all 103 sources and provide redshift constraints for 75 sources without 
previously known redshift. Out of those, eight have reliable photometric redshifts at z > 1 .3, while for the other 67 sources we pro- 
vide upper limits. Six of the former eight are BL Lac objects, which quadruples the sample of confirmed high-redshift BL Lac. This 
includes three sources with redshifts higher than the previous record for BL Lac, including CRATES J0402-2615, with the best-fit 
solution at z ~ 1.9. 

Key words. Techniques: photometric, (Galaxies:) BL Lacertae objects: general. Galaxies: distances and redshifts 


1. Introduction 

Since its launch in 2008, the Fermi Space Laboratory has dra- 
matically extended our view of the high-energy sky. The re- 
cently released 24-month catalog (2LAC) of Active Galactic 
Nuclei (AGN) detected by the Large Area Telescope (LAT; 
lAtwood et alJ 120091) revealed 885 high-sig nificance sources, 
the l arge majority of them being blazars (Ackermann et al. 
1201 lb . The latter form the most extreme class of AGN with 
their observational characteristics governed by the small an- 
gle between their relativisti c jets and the observer’s sight 
line felandfor d & Reeslll978l) . The resulting Doppler boosting 
makes blazars exceptionally bright sources at nearly all wave- 
lengths and therefore visible out to high redshift. 

The scientific relevance of blazars is very broad, ranging 
from laboratories for the physics and structure of relativistic jets 
(and thus the extraction of energy from the central massive black 


hole) (e.g.. lAbdo et al.|[2010dl) to probes of the extra-galactic 
background light (EBL) through attenuation of y-ray photons 
(e.g. . lAbdo et al.ll2010d) . One of the crucial parameters for these 
applications is the distance to the objects, which unfortunately is 
not easy to obtain in most cases. 

Two classes of blazars dominate the 2LAC population. These 
are the flat-spectrum radio quasars (FSRQs, 310 sources) and 
BL Lac objects (395), named after the prototype BL Lacertae 
dHoffmeisteill929l) . While for the former redshift measurements 
are routinely performed using their strong emission lines at UV- 
optical wavelengths, the featureless, power-law optical spec tra 
of BL Lac objects have proven to be a challenge (e.g JShaw et al.l 
120091) . Indeed, 220 of the 395 BL Lacs in the 2LAC (55 %) lack 
redshift estimates. Until this incompleteness is resolved, con- 
clusions about the EBL, the blazar sequence (e.g.. lFossati et al.l 
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19981 : Ghisellini et alii 19981) . and the blazar population in gen- 
eral remain tentative at best. 

Several methods have been exploited to increase the BL Lac 
redshift sample. At low distance, one can utilize the remarkably 
uniform absol ute brightness of t he giant elliptical BL Lac host 
galaxies (ISbarufatti et alJ 120051: iMeisner & RomanH 1201 Oh and 
very-high signal-to-noise optical spectroscopy to help to iden- 
tify weak emission or absorption features in a few other cases 
(e.g.. IShaw et al.l 12009 ). An alternative method, applicable to 
the more distant sources, is the photometric redshift technique, 
which consists of fitting spectral energy distribution (SED) tem- 
plates to multi-band photometry. Neutral hydrogen along the line 
of sight to the blazar will imprint a clear attenuation signature at 
the Lyman limit and thus allows an accurate estimate of the red- 
shift of the absorber. Even though the absorber will be located 
somewhere along the line of sight and its redshift will thus not 
necessarily correpond to that of the blazar, photometric redshifts, 
Zphot, will provide a reliable lower limit on the blazar redshift. 

In this paper we explore the use of ultra-violet to near- 
infrared quasi-simultaneous photometry to measure the redshift 
of 103 2LAC blazars, 86 of them without previous redshift 
constraints. The sample selection is described in Sect. [2] while 
the observations and data reduction are detailed in Sect. [3] and 
Sect. [4] The results and their discussion are presented in Sect. 0 
and Sect. 0 respectively. 

2. Sample Selection 

The sample is selected mainly from the clean sample 
of the 24-month catal og of AGN detected by Fermi/ LAT 
dAckermann et al.ll201 lh which includes 885 high-significance 
y-ray sources which are statistically associated with AGNs and 
located at high Galactic latitudes ( \b\ > 10 deg). To accomodate 
the location of one of our follow-up instruments in the southern 
hemisphere, we applied a cut in declination ( t ) .12000 < +25 deg). 
Furthermore, the importance of the ultra-violet bands for the 
Zphot estimates lead us to discard all sources with Galactic fore- 
ground reddening of E b-v > 0.2 mag, as derived from the 
ISchlegel et al.l ( 1 998 ) maps. The final sample is composed of 80 
sources without known redshift but with optical or radio counter- 
part associations from Bayesian statistics (lAbdo et al .1(20 1 O alihl). 
Likel ihood Ratio, or log N - log S methods ( Ackerman n et al.l 
1201 lh and was further extended by eight sources that do not be- 
long to the clean 2LAC sampkfl In addition, observations for 16 
2LAC sources with existing redshift measurements are included 
for verification. We emphasize here that our sample is not sta- 
tistically complete, but biased by selection. The final list of 104 
targets is presented in Table Q] Of those, 82 have been classified 
as BL Lac obje cts, three as a FSRQ, an d the remaining 19 are of 
unknown type dAckermann et alJl201 ll) . 



0.01 0.10 1.00 10.00 100.00 

^grond-uvot [days] 


Fig. 1. Cumulative distribution of the (absolute) time between 
the mid-points of the GROND and Swift/ UVOT observations. 
The median offset is 1.5 d. For longer UVOT observations, 
which were spread over several satellite orbits, the mid-time of 
the nvw2-exposure was taken as reference. 

The UVOT pointings have been obtained as part of the Swift 
fill-in program from January 2010 to October 2011. With few 
exceptions, they comprise observations in three optical ( u , b, v) 
and three UV (uvw2, uvm2, uvw\ ) lenticular filters, covering the 
wavelength range between 160 nm and 600 nm (iPoole et al.l 
120081) . The exposure times vary from source to source de- 
pending on visibility and brightness of the optical counterpart. 
Approximate values are 100 s each in u, b, v and 200, 250, 400 s 
in uvwl,uvm2,uvw2, respectively. For a few sources multiple 
consecutive orbits, and thus longer exposures, were obtained. 

While UVOT needs to cycle through all filters in sequence, 
GROND observes in it’s four optical g ' , r ' , i ' , 7 / and three near- 
IR J, H, K s bands simultaneously. This capability is particularly 
important for studies of fast-varying sources such as blazars as it 
allows the construction of a reliable spectral energy distribution 
(SED) without the need to correct for source-intrinsic variabil- 
ity. Unfortunately, scheduling limitations rarely allowed simul- 
taneous observations with both instruments, UVOT from space 
and GROND from the ground. However, in order to minimize 
the impact of variability on the combined SED, GROND obser- 
vations were executed as close as possible in time to the Swift 
pointings. Except in a few cases where visibility or weather con- 
straints lead to offsets of more than 10 d, observations were typ- 
ically achieved within 1-2 d of each other (see Fig. The Zphot 
uncertainties associated with these offsets will be discussed in 
Sect. 1431 

A typical GROND observation had an integration time of 
2.4 min in (g', r' , i',z') and 4.0 min in ( J,H , K s ). 


3. Observations 


The observations were performed with the Ultraviolet and 
Optical Telesco pe (UVOT; Roming et al.l 120051) onboard the 
Swift satellite dGehrels et alJ l2004li and the Gam ma-ra' 
Optical/Near-Infrared Detector (GROND; iGreiner et alJ 1200 
mounted on the MPI/ESO 2.2 m telescope at La Silla, Chile. 
Each source has been observed at least once with both instru- 
ments, although for some targets several epochs had been ac- 
quired. 


4. Data & Analysis 

4. 1 . Swift/UVOT and GROND 

UVOT photometry was carried out on pipeline processed sky 
images downloaded from the Swift data center Q, following the 
standard UVOT procedure dPooie et al.l 12008 ). Source photo- 
metric measurements were extracted from the UVOT imaging 
data using the tool uvotmaghist (vl.l) with a circular source ex- 
traction region that ranged from 3"5-5" radius to maximise the 


1 flagged due to analysis issues by the LAT team 


2 http://www.swift.ac.uk/swift_portal 
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Table 1. Target List (first five rows, rest in Journal only) 


Name 

FGL Name” 

0112000* 

A b 

0 J2000 

offset ["] c 

E B -v [mag]" 

CRATES J000 1-0746 

2FGLJ0000.9-0748 

00:01:18.00 

-07:46:26.9 

0.2 

0.03 

CRATES J0009+0628 

2FGLJ0009. 0+0632 

00:09:03.93 

06:28:21.3 

1.1 

0.07 

CRATES J002 1-2550 

2FGLJ002 1.6-2551 

00:21:32.56 

-25:50:49.1 

0.8 

0.02 

1RXS J002209.2- 185333 

2FGLJ0022.2- 1853 

00:22:09.26 

-18:53:34.9 

3.1 

0.03 

RX J0035.2+1515 

2FGLJ0035.2+ 1515 

00:35:14.71 

15:15:04.2 

3.4 

0.07 


a : FennifLAT ID. 

b : GROND coordinates of the optical counterpart with typical uncertainties of (y.'3 in both directions. 
c : Offse t of counterpart locati on from those given in Ackerm ann et all d2011h . 
d : FromUchlegel et al] d 1 9981) . 


signal-to-noise. In order to remain compatible with the effec- 
tive area_calibrations, which are based on 5” aperture photome- 
try d Poole et al.l2008 1. an aperture correction was applied where 
necessary. This correction was at maximum 5-6 % of the flux, 
depending on the filter. For the further analysis, all magnitudes 
were converted into the AB system and are presented in Tab. [2] 
Typically achieved 3-cr limiting magnitudes are uvw2 AB ~ 22.4 
and b AB « 19.6. 


The GROND data (Tab. 0 were reduced and analyzed 
with the standard tools and methods described in iKriihler et alJ 
d2008l) . The g'.r',i',z' photometry was obtained using point- 
spread-function (PSF) fitting and calibrated against observations 
of fields covered by the SDSS Data Release 8 (iAihara et al.l 
1201 ll) . Due to the undersampled PSF in the near-infrared, 
the J, H, K s photometry was measured from apertures with 
the sizes corresponding to the Full- Width at Half Maximum 
(FWHM) of field stars and calibrated against selected 2MASS 
stars (tSkrutskie et al.l [ 200 6). This resulted in lcr accuracies of 
0.04 mag (g',z'), 0.03 mag (/, i'), 0.05 mag ( J,H ), and 0.07 mag 
(K s ) for the absolute calibration. While the SDSS calibration di- 
rectly provides magnitudes in the AB system, the near-IR pho- 
tometry required further transformation from the 2MASS -native 
Vega system into AB. Typically achieved 3-cr limiting magni- 
tudes are r' AB 23.5 and Ws.ab ~ 19.8. 

Correction for Galactic foreground extinction wa s perfo rmed 
following the procedure describ ed in lCardelli et al.l (fl9891 with 
E b-v from lSchlegel et akl (11998 ). For the UVOT bands, the cor- 
rection factors presented in iKataoka et al. (120081) were used. 
Uncertainties in E b-v (10%; [Schlegel et al.ll998 ) and in the red- 
dening law produce an additional contribution to the photometric 
error budget. As this systematic uncertainty is coupled between 
the the photometric bands for a given position in the sky, its im- 
pact on the SED fitting is smaller than the contribution to the 
photometry in each individual filter. An exact calculation is very 
complex. Thus, we adopt a conservative 5 % of the reddening 
value in each band which is added in quadrature to the photo- 
metric uncertainties. 


4.2. Counterpart selection and morphology 

Astrometric solutions for the GROND optical bands have been 
obtained through comparison with USNO or, if available, with 
SDSS measurements, achieving a typical rms of 073 in both co- 
ordinates. For the large majority of our targets (91) only a single 
optical source was found within 2" radius of the location given 
in the 2LAC catalog. For an additional twelve objects isolated 
candidate counterparts were detected at offsets of up to 6". In 
one case (CGRaBS 1407-4302), several blended sources made 
a reliable identification impossible. Counterpart coordinates and 



r’-band Brightness [mag] 


Fig. 2. Observed GROND /-band magnitude distribution of the 
optical counterparts. 


offsets for the 103 sources with identification are given in Tab.0 
Most of the identified optical sources are bright and point-like 
with typical magnitudes in the GROND /-band of 16-18 mag 
(Figure 0. 

The source morphology can be used as a prior for the se- 
lection of the SED template library (see Sect. 14.41) and redshift 
range, i.e., clearly extended sources will have a contribution 
from the host galaxy emission to the observed SED and are likely 
located at low redshift. A simple morphological distinction in 
point-like and extended has been performed on the GROND im- 
ages and will be used for discussing the reliability of the z ph0 i 
estimates. 


4.3. Variability correction 


The optical emission of blazars is known to vary on a wide 
range of time scales. Brightness variations range from_a few 
tenths of a magnitude within minutes to days (e.g. 


fl970L iMiller et all [[989; Ca rini et all 179911 llJrrv et al 


Racim 


19931 


to several magnitudes ov er weeks to years (e.g., Cinrini et al.l 
120031: iRaiteri et al.l 120051) . This variability contributes signifi- 
cantly to the uncertainties when constructing an SED from non- 
simultaneous multi-band, multi-instrument observations. 

As GROND observes in all seven bands at the same time, the 
g' to K s photometry can be considered as a snapshot of the SED 
and is thus unaffected by variability on time scales longer than 
the exposure of an individual observation. This leaves two ar- 
eas where a proper treatment of variability has to be performed, 
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Table 2. Swift/UV OT photometry (first five rows, rest in Journal only) 


Name 

UT Date" 



AB Magnitude 1 ’ 



A mcR->uv c 



uvwl 

uvml 

uvwl 

u 

b 

V 

[mag] 

CRATES J0001... 

2010/10/14 10:40 

18.99 ± 0.08 

18.74 + 0.10 

18.58 + 0.09 

18.18 ±0.09 

17.88 + 0.10 

17.90 ±0.19 

0.01 

CRATES J0009... 

2010/12/07 08:55 

20.16 + 0.15 

19.71 ±0.19 

19.69 + 0.18 

19.24 ±0.19 

18.51+0.18 

18.14 ±0.24 

-0.24 

CRATES J0021... 

2010/11/18 18:32 

18.86 + 0.08 

18.90 + 0.27 

18.57 + 0.09 

18.15 ±0.10 

17.82 + 0.11 

17.93 ±0.21 

0.09 

1RXS J002209.2... 

2010/11/18 13:49 

18.36 ± 0.07 

17.95+0.08 

17.91+0.08 

17.48 ± 0.07 

17.20 + 0.08 

16.97 ±0.12 

-0.08 

RX J0035.2... 

2011/01/14 01:08 

18.36 + 0.08 

17.95+0.08 

17.91+0.08 

17.48 ± 0.07 

17.20 + 0.08 

16.97 ±0.12 

-0.24 

Approximate start time of uvwl exposure. 








. Corrrected for Galactic foreground reddening. Upper limits are 3-cr. 






c : Variability-correction factor to be applied to GROND photometry (see text). 





Table 3. GROND photometry (first five 

rows, rest in Journal only) 






Name 

UT Date" 




AB Magnitude 6 






S' 

r' 

i ' 

z 7 

J 

H 

K s 

CRATES J0001... 

2010/10/14 03:33 

17.80 + 0.05 

17.43 ± 0.05 

17.18 + 0.05 

16.86 ± 0.05 

16.48 ± 0.05 

16.16 ±0.06 

15.76 + 0.08 

CRATES J0009... 

2010/12/12 00:41 

18.91 ±0.05 

18.49 ± 0.05 

18.15 + 0.05 

17.89 ± 0.06 

17.34 ± 0.06 

16.90 ± 0.06 

16.46 ± 0.08 

CRATES J0021... 

2010/11/23 00:13 

17.69 ± 0.05 

17.43+0.05 

17.20 + 0.05 

16.98 ± 0.05 

16.76 ± 0.05 

16.51 ±0.05 

16.15 + 0.07 

1RXS J002209.2... 

2010/11/23 00:21 

17.30 ± 0.05 

16.92 ± 0.05 

16.71+0.05 

16.51 ±0.05 

16.22 ± 0.05 

16.00 ± 0.06 

15.57 ± 0.08 

RX J0035.2... 

2011/01/14 00:48 

17.39 ± 0.05 

17.10 + 0.05 

16.96 + 0.05 

16.77 ± 0.05 

16.41 ±0.06 

16.17 ±0.05 

15.92 ± 0.08 


Exposure start time. 

b \ Corrrected for Galactic foreground reddening. Not corrected for variability. Upper limits are 3-cr. 


i) between the exposures in the individual UVOT filters and ii) 
between the GROND and UVOT pointings. 

UVOT observes in each band separately and, for our pro- 
gram, typically cycled through its filters in a specific order 
( uvwl , u, b, uvw2, v, mm2). A complete sequence in all six bands 
took « 12 min, during which time the target may vary in bright- 
ness. Without a priori knowledge of the spectral shape, this vari- 
ability can only be accounted for statistically, i.e., by including 
an additional contribution to the systematic uncertainties. For 
this purpose we analyzed the photometry for those blazars for 
which multiple UVOT exposures in a single band were available. 
Here, we focused on the bluest filter (uvwl) which would be 
least affected by a potential host galaxy contribution and should 
thus give a good representation of the maximum variability. In 
Fig. 0 the distribution of variability as function of time between 
two uvw2 exposures is shown for 25 targets with accurate pho- 
tometry (6m uvw 2 <0.1 mag) in each image. While variability as 
large as Am UVW 2 ~ 0.4 mag can be seen for some sources also on 
time scales of At < 0.1 d, the median variation for At < 50 d is 
« 0.1 mag, comparable to the photometric uncertainty. 

The shortest time sampled between two nvw2 exposures in 
this sample is one satellite orbit (~ 90 min) and therefore longer 
than it takes to cycle through the six filters. However, it is rea- 
sonable to assume that the median variation on a 12 min time 
scale is not much larger than the one over 96 min. Thus, we con- 
servatively adopted Am =0.1 mag as an additional contribution 
to the systematic uncertainties for all UVOT filters. 

In order to correct between GROND and UVOT we made 
use of the_spectral overlap provided by the two instruments 
(iKriihler et al] 1201 ll) . Under the assumption that the SED re- 
mains unchanged and can be approximated by a power law and 
thus, that the GROND and b, v photometry follow the same spec- 
tral slope Q color-terms can be used to derive the normaliza- 


3 This is reasonable for redshift z < 3.5 and for sources where the 
host galaxy contribution to the emission is negligible. See Sect. l4.4l for 
further discussion. 



At [days] 

Fig. 3. Absolute brightness changes in the uvw2 filter as func- 
tion of time between two UVOT pointings for all 25 sources 
with multiple observations. The black crosses show the median 
values (and median absolute deviation) and the grey dots mark 
all 174 pairs of observations. While micro-variability (A m UVW 2 ~ 
0. 1 mag) dominates on time scales shorter than 50 d, signif- 
icantly larger variability (Am UVW 2 > 0.3 mag) is observed at 
Af > 50 d. 

tion offset. Here we used the calibration with the largest spectral 
overlap, namely the one between b,g' , and r'\ 

b-g' =Q.\l{g' -r') + 0m{g' -r'f (1) 

which is valid for -1 < ( g ' - r') <2 (IKriihler et al.ll201 1). 
These corrections, AmagGR->uv (Tab [2}, were calculated for all 
sources with accurate photometry in b, g\ and r' (6m <0.2 mag). 
They were computed for each object individually and applied to 
the GROND photometry, shifting it in line with UVOT. Typically 
values are of the order of +0.1 mag (Fig. 0J. In cases where 
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At [days] 

Fig. 4. Absolute brightness offset as function of time between the 
GROND and Swift/UV OT observations. The black crosses show 
the median values (and median absolute deviation) and the grey 
dots mark all 90 sources for which the offset was determined 
(excluding eleven sources with large photometric uncertainties, 
missing b, g', or r' photometry, or strong host galaxy contribu- 
tion). 


multiple epochs were available, we used the closest pair of 
GROND and UVOT observations. The final SED for each source 
corresponds to the best possible reconstruction between 160- 
2200 nm at the time of the UVOT observation. No correction 
was performed for sources where the photometric uncertainty 
(6m > 0.2 mag in b,g', and/or r'j was larger than the typical 
Ama.gGR-.uv- 

For a few sources, the emission from the underlying 
host galaxy contributes significantly to the overall SED (see 
Sect. 14.41) . While the blazar emission can still vary, the host con- 
tribution naturally remains unchanged, violating the assumption 
of a constant SED slope. Thus, no variability correction was per- 
formed in those cases where the UV to optical SED clearly de- 
viates from a power law. 

4.4. SED fitting 

The photometric redshift computation was performed with the 
publically availa ble LePhare code v.2.2 (iArnouts et al] fl999t 
lllbert et al.l2006l) . The program uses the simple^ 2 fitting method 
to differentiate between a theoretical and an observed photo- 
metric catalogue. We used three custom template libraries for 
comparison with our GROND and UVOT photometry. The first 
was composed of 40 power law SEDs of the form F.i oc 
with fi ranging from 0 to 2 in steps of 0.05 and describes pri- 
marily the central-engine dominated blazars. Here, the source- 
intrinsic extinction was assumed to be zero and the luminos- 
ity prior was choosen to cover the expected range of absolute 
magnitudes (-20 < Mg < -30: I Veron-Cettv & Veronll2010k . 
Such a simple model is a good approximation for the UV-near- 
IR SEDs of BL Lacs, which form the bulk of our sample. Also 
FSRQs, if their emission lines are not dominant, can be mod- 
eled with a power law to first order. The second library con- 
tains tem plates of normal (inactiv e) galaxies and galaxy/AGN 
hybrids (ISalvato et al.ll2009l.l201 lh and has been implemented in 
order to model host galaxy dominated objects. If the host con- 
tribution is significant, the 4000 A break can emerge and poten- 
tially be misinterpreted as a Lyman-limit. While for the galaxy 


models intrinsic extinction was again neglected, a different lu- 
minosity prior of -8 < Mg < -30 was chosen to allow for the 
expected fainter host magnitud es. The third library consists of a 
wide range of stellar templates dPickles 1998; Bo hlin et al.l l995: 
Chab rier et al]|2000h . included to test against potential false as- 
sociations. For each source in our sample, all three libraries were 
fit independently. 

The most prominent spectral feature used to measure z p h 0 t 
for blazars is caused by the absorption through neutral hydrogen 
along the line of sight. At z ~ 0.8 the Lyman-limit moves into 
the wavelength covered by the UVOT and starts to suppress the 
flux in the uvw2 filter. While this is the theoretical lower limit on 
Zphot for the combined GROND and UVOT coverage, uncertain- 
ties associated with the photometry and modeling will shift it to 
a higher redshift. Usually one would use a spectroscopic train- 
ing sample to evaluate the redshift range and accuracy which is 
accessible with our SEDs. As the majority of the sources with 
known redshift in our sample are located at z < 1, they alone 
do not suffice to provide a statistically meaningful test of the 
method. Instead we assess the reliability in our photometric red- 
shift computation with a Monte-Carlo approach. 

Here, we simulated 27,000 test power-law SEDs with spec- 
tral slopes ranging from 0.5 to 2, redshifts from 0 to 4, and at 
three apparent magnitudes, r' = 17, 18, 19mag. For each source, 
the photometry in the individual bands was scattered around the 
model magnitudes with a brightness-dependent statistical con- 
tribution and with a component representing the systematic un- 
certainties from the photometric calibration and variability cor- 
rection. The resulting SEDs were fed back into LePhare and the 
output Zphot was compared to the input value. 

As can be seen in Fig. El there is a good correspondence be- 
tween the input and the recovered redshift for z s i m > 1.2. In par- 
ticular for bright sources (/ < 18 mag) the number of outliers rj, 
defined as the number of sources with |Az(l +z s i m )l > 0.15, drops 
steeply above this redshift. For fainter sources (r' ~ 19 mag), the 
increasingly less constrained UV photometry causes the outlier 
fraction to remain high (~ 20%) out to z S j m ~ 2. A more quan- 
titative selection can be made when using P z , the integral of the 
probability distribution function f f(z)dz at z p h 0 t ±0.1(1+ z p hot), 
which describes the probability that the redshift of a source is 
within 0.1(1 + z) of the best fit value. When choosing a cut at 
P z > 90%, the majority of the outliers at z p h 0 t > 1-2, as well as 
nearly all the solutions with z p h 0 t <1.2 disappear (inset Fig. [5}. 
This confims that for power-law sources with z < 1 .2 the avail- 
able photometry can only place an upper limit, and suggests that 
P z >90% can be used as a criterion to identify reliable photo- 
metric redshift solutions. Selecting a higher P z threshold, e.g., 
99 %, can further reduce the outlier fraction but will also shift 
the lower constrainable redshift bound to higher values. In order 
to reach such a very tightly peaked redshift probability distribu- 
tion function more than one photometric band has to be affected 
by the Lyman-limit, limiting the redshift range to z > 1.5. 

5. Results 

The fit results for all 103 sources with identified counterparts 
are presented in Tab. [4] Here we give the photometric redshifts 
(and their 90 % confidence errors), the P z values, and the best 
fit models for the power law and galaxy templates. No source 
requires a stellar template. 

As discussed in the previous section, P z > 90 % can be used 
as a good reliability criterion for the photometric redshift solu- 
tion. Applying this cut to the power law model fits results in 
a sample of 15 sources with x 2 < 30. Next, we also require 
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Table 4. Excerpt of the results table. Only entries for those sources with reliable redshifts of z p hot >1.2 are shown, (rest in Journal 
only) 


Name 


£phot,best 


1-28!";" 
1 44+0-I6 

1 / —0.17 

I 47+o.n 

92^ 

' J 

a ' jj -0.13 

1 54 +014 
— Q-19 

1 60 +01 ° 
1 - u -0.05 

i 77+0.13 

rifiM 

-m 
0.18 



£phot 

0 79+ ul4 
Q4j+8:J1 
-QQ8 

0 47 + 0-23 
' -0.47 

l.ll+o-io 

L28 -qI 

1-58!™ 

0 46 +0 ^ 4 

u ‘^-0.05 

1 i q+0.50 

j 54+8:18 

-Q.Q9 

1 60 + °' 4 
-QQ9 

0 55 + "-' 6 

-0-14 


X 


RX J0035.2+1515 
PKS 0047+023 
PKS 0055-328 
PKS 0332-403 
CRATES J0402-2615 
SUMSS J053748 571828 
PKS 0600-749 
CRATES J0630-2406 
OM 235 

CRATES J 13 12-2 156 
CLASS J2352+ 1749 


1.45- 


1.426 


1.549 

1.491 


8.7 

13.6 

11.2 

13.0 

9.5 

21.5 

7.5 
24.6 
11.9 
18.4 
7.1 


a : Best photometric redshift, 
b : Spectroscopic or imaging 
c : Photometric redshift with 90% confidence uncertainties. 


see text. 

redshift (if known) from l Ackermann et all d201 li) . 


galaxy 

P z c model 

"751 12249 1.60.TQSO 1.40 

80.5 pLQSO_DR2_029_tO 

38.4 12249 1.50.TQSO 1.50 

62.3 12249 1.60.TQSO 1.40 

48.9 12249 1.60.TQSO 1.40 

99.6 pl.I22491_30.TQS01_70 

58.1 12249 1.60.TQSO 1.40 

83.7 12249 1.60.TQSO 1.40 

94.9 pl.I22491_30.TQS01_70 

100.0 12249 1.60.TQSO 1.40 

64.2 pl_I22491-20_TQSO1.80 


d : Redshift probability density at z phot ±0.1(1 + z phot ). 
e : Spectral slope for power law model of the form 

X. Starburst/QSO hybrid templates with varying contributions of the two components. Template names reflect the relative contribution of the 
starburst (1 2249 1.NN) and AGN (TQSOIJMN) models to the hybrid, “pi” marks an additional power-law component at short wavelengths. See 
ISalvato et all J20 1 ll) for details on the templates. 


that a source can not be fit well with a galaxy /AGN-hybrid tem- 
plate at low redshift ( z < 1.2). This step is required in order 
to identify those sources for which a degeneracy in the pho- 
tometric redshift solution does not allow a distinction between 
a low-redshift galaxy/AGN-hybrid and a blazar at higher red- 
shift. In analogy with our selection for a good redshift constraint 
(P z > 90 %), we consider P z , ga i < 90 % as a reliable criterion 
that no adequate low-redshift galaxy/AGN-hybrid solution was 
found. This removes further seven sources from the sample and 
leaves eight candidates, all unresolved in the GROND images 
and with z pho t, P i > 1-2. 

For three additional sources the power law and galaxy mod- 
els fit similarly well (P z , p i ~ P Zlga i > 90 %) and give compara- 
ble high-redshift solutions, indicating that the Lyman-limit was 
fit by power law and galaxy templates alike. Including those, 
we end up with a sample of eleven sources with reliable photo- 
metric redshift computation, all , with - ph()L ranging from z phol ~ 
1.28 to z p hot ~ 1.92. This encompasses eight sources without 
previous redshift measurements and three blazars with known 
redshift. For the latter, the photometric redshifts are within 3- 
cr of the known value (see Fig. O, indicating an accuracy of 
Az/(1 + Zs P ec) < 0.12 for our z p h 0 t method described in Sect. 14.41 
The fit results and SEDs for the eleven sources are presented in 
Tab.[4]and Fig. [6] respectively. 

For 21 sources, the galaxy library provides the better fit with 
Pz.gai > 90%, P ZiP i < 90%, and x\ & \ < 30. All of these are 
best described by starburst/QSO hybrid templates at z < 1-3 
and none require an elliptical galaxy, i.e., a strong 4000 A break. 
The most significant difference between the starburst/QSO hy- 
brid templates and the power-law models is that the former ex- 
hibit^ strong e missi on lines originating from star-formation (see 
ISalvato et al.l I2009L for more details on the templates). These 
lines can be fit to small deviations in the photometry from a 
power-law and thus provide low x 2 values and apparently well- 
constrained z p hot solutions. However, as described in Sect. 14.31 
the variability correction is not reliable when applied to non- 
power-law sources, as the color terms, and thus the correction 
factors, will be erroneous. Also, 18 of these sources have pre- 


viously been classified as BL Lac objects and are therefore 
unlikely to show strong emission lines in their optical spec- 
tra. Given these considerations, we do not regard their galaxy- 
template photometric redshifts to be reliable. 

There are 79 sources which are well-fit by power-law tem- 
plates (x 1 < 30) but have a low P z . In other words, the spec- 
tral slope is well defined but the photometric redshift cannot 
be unambiguously determined. As expected, the best solutions 
are predominantly at low redshift (z < 1.2) with typical red- 
shift probability distributions that are flat down to z = 0. For 
those we give the 90% upper limit of the photometric redshift.. 
This sub-sample also includes eleven sources with a previously 
known redshift from spectroscopy or imaging, all in agreement 
with our limits. 

We also provide the 90 % upper limit for four sources for 
which the power-law and galaxy libraries give similarly good 
constraints (P ZiP i »P z , ga i > 90%) but with very different red- 
shift solutions. This degeneracy can occur, e.g., when a break 
in the SED is similarly well fit by the Lyman limit of a power- 
law model and with the 4000 A break of a galaxy template. No 
Z P hot,best is given for ten sources for which no satisfactory fit 
( x 2 > 30) was obtained. 


6. Discussion & Conclusion 

In this paper we presented redshift constraints for 103 blazars 
from the 2LAC catalog using UV-to-near-IR multi-band pho- 
tometry obtained quasi-simultaneously with Swift fUWOT and 
GROND. We provided the first reliable redshift measurements 
for eight sources and new upper limits for an additional 66 tar- 
gets. Of the eight sources with reliable redshift, seven are lo- 
cated at z p hot > 1-3. Six of those belong to the BL Lac popula- 
tion. For comparison, out of the total 395 BL Lac in the 2LAC 
sample only two sources have previously been known to lie at 
Z > 1.3 (Fig, [7}. Redshifts for these two, PKS 0332-403 and 
CRATES 113 12-2156, have also been confirmed with our pho- 
tometric observations. 
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Fig. 5. Recovered best-fit photometric vs input redshift (top), ac- 
curacy (middle), and outlier fraction (bottom) for 27,000 sim- 
ulated sources with spectral slope of 0.05 < /3 < 2 and op- 
tical brightness of r' = 17, 18, 19mag. The pink circles indi- 
cate the location of three blazars in our sam ple with known 
redshift, PKS 0332-403 (z spec = 1 .35 1 . [Bergeron et all l20H . 
OM 2 35 (tw - 1.549), and CRATES J1 3 12-2156 (z spec = 
1 .49 1 . 1 Ackermann et al. 20111), as well as that of PKS 0537- 
286 (z spe c = 3.104, Wright et al.l 1 19781) . a LAT-non-detected 
blazar with available GROND and UVOT photometry. The green 
area in the middle panel corresponds to the forbidden parameter 
space where z p h 0 t > 0. Systematics are visible at z p h 0t < 0.8. 
The large majority of the fits have < 30 (for typically 10 
d.o.f.) with the exceptions shown in light grey. The inset in the 
top panel displays the results if only solutions with P z > 90 % 
are selected. 


The six new BL Lac redshifts at z > 1.3 represent a dra- 
matic (from two to eight) increase of the confirmed high-z Fermi 
sample and thus demonstrate the opportunity that the SED tem- 
plate fitting technique holds for obtaining photometric redshifts 
for BL Lac sources. Lor our sample, this was possible due to 
the densely-covered wide spectral range (1 60-2200 nm), neces- 
sary for a reliable constraint of the spectral slope, the excellent 
ultra-violet coverage to measure the Lyman-limit, and the quasi- 
simultaneity of the observations, important for minimizing the 
impact of source-intrinsic variability. The method applied in this 
work overcomes significant challenges inherent in other redshift 
techniques, namely the simplicity of the optical emission of BL 
Lacs manifested as the power-law-shaped synchrotron spectral 
component. This complicates spectroscopic redshift measure- 
ments, which are in most cases limited to the optical wavelength 
regime and thus insensitive to the Lyman-limit. Instead, spec- 



Fig. 6. SEDs for the eleven sources for which photometric red- 
shifts larger than 1 .2 were estimated (see Tab. 0} and one ex- 
ample for z p hot < 1-2. From bottom to top: CRATES J0402- 
2615, CRATES J1312-2156, OM 235, CRATES J0630- 
2406, SUMSS J053748-571828, PKS 0600-749, PKS 0332- 
403, CLASS J2352+1749, PKS 0047+023, PKS 0055-328, 
RX J0035. 2+1515, and BZB J0543-5532. 


troscopy relies on the detection of very faint emission features, 
mainly from the underlying host galaxy. It therefore requires 
very high signal-to-noise, often at the cost of long exposure 
times. On the other hand, the photometric redshift method for 


14 - 



0.5 1.0 1.5 2.0 

z 

Fig. 7. Redshift distribution for the nine BL Lac objects with re- 
liable photometric redshifts (blue filled histogram) together with 
the dist ribution for the BL Lac in the 2LAC catalog (empty his- 
togram ;jAckermann|eran20Tl|). Sources with z < 0.5 have been 
omitted for clarity. 
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BL Lac does not require significant telescope time, in particular 
when data from efficient multi-channel instruments like GROND 
are available. Its weakness, however, is that the method relies 
on the detection of a particular spectral feature, which, albeit 
strong, is located in the rest-frame far-ultra-violet. For z, < 3, this 
is only accessible with space-based ultra-violet telescopes, and 
even then usually limited to redshifts above z ~ 1.2, as shown 
above. 

The sensitivity of this method to high-redshift sources also 
allows the placement of upper limits for those SEDs that do not 
show the imprint of the Lyman-limit. Only for three targets is 
this boundary above z = 2. In one case, CRATES J0250+1708 
(z < 3.1), the optical counterpart is too faint (r' » 20.7 mag) 
and only upper limits in all six UVOT bands could be obtained. 
The counterpart of CRATES J0705-4847 is only detected in the 
GROND J (« 19.7 magAB) and H bands and its SED and redshift 
are thus poorly constrained. Finally, the SED of ATG20 J0124- 
0625 (z < 2.46) shows a significant break at z ~ 1.9 + 0.5. 
However, as the redshift probability distribution is broad, and 
the x 2 of the power law fit comparable to that of a z ~ 0.4 so- 
lution for a galaxy template, the photometric redshift is consid- 
ered to be unreliable. Except for those three, no other source in 
our sample has a best-fit photometric redshift at z > 2. Thus, 
CRATES J0402-2615 can now be considered the most distant 
known BL Lac with a measured redshift of z ~ 1 .92. 

Fig Hindicates that our high-redshift findings are the natural 
extension of the existing 2LAC redshift sample. However, due to 
the incompletenesses of both samples, we refrain from drawing 
any quantitative conclusions at this stage. We note, however, the 
our resul t is g eneral agree ment with the theoretical predictions 
from lCiommi et al. ( 201 1). A more detailed physical interpreta- 
tion of the establishment of an increased fraction of high-z BL 
Lacs in the 2LAC sample will be reported separately. 
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